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ABSTRACT 

We have studied the dynamical evolution of rotating star clusters with mass spectrum 
using a Fokker-Planck code. As a simplest multi-mass model, we first investigated 
the two-component clusters. Rotation is found to accelerate the dynamical evolution 
through the transfer of angular momentum outward, as well as from the high masses to 
the low masses. However, the degree of acceleration depends sensitively on the assumed 
initial mass function since dynamical friction, which generates mass segregation, also 
tends to accelerate the evolution, and the combined effect of both is not linear or mul- 
tiplicative. As long as dynamical friction dominates in the competition with angular 
momentum exchange the heavy masses lose random energy and angular momentum, 
sink towards the centre, but their remaining angular momentum is sufficient to speed 
them up rotationally. This is gravo-gyro instability. As a consequence, we find that the 
high mass stars in the central parts rotate faster than low mass stars. This leads to 
the suppression of mass segregation compared to the non-rotating clusters. From the 
study of multi-component models, we observe similar trends to the two-component 
models in almost all aspects. The mass function changes less drastically for clusters 
with rotation. Unlike non-rotating clusters, the mass function depends on R and z. 
Our models are the only ones that can predict mass function and other quantities to 
be compared with new observations. 
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1 INTRODUCTION 

Since the pioneering studies of Lupton & Gunn (1987) and 
its subsequent applications to fit actual globular cluster ob- 
servations it is evident that galactic globular clusters exliibit 
some degree of rotation, and tliat tliere is also a consis- 
tent amount of observed flattening of their shapes (White & 
Shawl 1987). More recently, van Leeuwen et al. (2000) and 
Anderson & King (2003) measure globular cluster rotation 
in proper motion, in order to derive the true direction of its 
rotational axis. Though not dominant the amount of rota- 
tional energy could also not be neglected on the other hand 
side, as was seen by the models of Einsel & Spurzem, which 
showed that even a moderate fraction of rotational energy 
in a cluster leads to significantly faster core evolution. 

This is the third in a series of studies on the dynamical 
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evolution of rotating stellar system by using orbit-averaged 
2D Fokker-Planck (FP) models that include the effects of 
initial rotation. In the previous two papers, pre-coUapse evo- 
lution (Einsel & Spurzem 1999, Paper I) and the evolution 
after core-collapse (Kim et al. 2002, Paper II) of rotating 
clusters composed of equal mass stars were studied. In the 
present study we explore the dynamical evolution of the ro- 
tating stellar systems with mass spectrum as a natural ex- 
tension to the previous works. We expect that the exchange 
of angular momentum between different mass species would 
significantly affect the course of dynamical evolution, as the 
energy exchange is known to play important role in the evo- 
lution of multi-mass clusters. 

Direct integration of the Fokker-Planck equation is used 
as a statistical method. Comparisons between results ob- 
tained with FP method and results from A'^-body simulation 
show that the approximations and assumptions which were 
used in FP models are reasonable, but need to be checked 
carefully by comparison of different methods. The simplest 
implementation of an FP model is a one-dimensional (ID) 
FP model where the distribution function / is assumed to 
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depend only on the energy E of stars. All physical prop- 
erties depend on the distance from the cluster center only 
(Cohn 1980; Lee, Fahlman & Richer 1991). The use of ID 
models is inspired by the fact that the shape of globular clus- 
ters is approximately spherical, but it ignores the anisotropy 
of the velocity dispersion. Observations of globular clusters 
and theoretical models suggested that there exists a differ- 
ence in the velocity dispersion between radial and tangen- 
tial direction, especially for stars in the outer regions of the 
clusters (Takahashi 1995, 1996, 1997). The two-dimensional 
FP model where distribution function depends on energy 
and angular momentum with spherical geometry has been 
pioneered by Cohn (1979) but studied extensively only re- 
cently because of difficulties in rmmerical integration of the 
equations (Takahashi, Lee & Inagaki 1997, Takahashi & Lee 
2001). In our case, modelling axisymmetric systems with 
only two integrals requires a careful check in particular. Any 
third integral is neglected here completely, which means that 
the diffusion properties of orbits in the axisymmetric poten- 
tial are treated as a function of two integrals only. 

Although there are successes in previous FP models to 
explain the dynamical evolution of star clusters, only few of 
them considered the natural and important physical prop- 
erty of the existence of an initial angular momentum in the 
star cluster. Kim et al. (2002) have reviewed briefly the pre- 
vious studies where the initial rotation is considered. They 
presented the first post-core collapse studies on the evolution 
of rotating stellar systems and found that the global shape 
of the rotational structure of the globular cluster changes 
little, though the strength of the rotation (measured using 
the magnitude of the rotational velocity or the «-componcnt 
of the angular momentum) decreases continuously with time 
due to the outward transfer of the angular momentum. An 
incorporation of mass loss and enhanced two-body relax- 
ation processes accelerated the evolution of the star cluster 
not only in corc-coUapso time, but also the dissolution time 
of the cluster. They found an approach to self-similar evo- 
lution in late core- collapse. 

The previous studies where the initial rotation is consid- 
ered, assumed a star cluster with equal mass stars. Inclusion 
of the initial mass spectrum known to change the evolu- 
tion of star clusters significantly, shortening of core-collapse 
times for example (Lee, Fahlman & Richer 1991). We also 
expect that there may be important physical processes be- 
tween different mass species concerning the exchanges of the 
angular momentum. Until now, no study has been done on 
rotating star clusters with mass spectrum except for one pre- 
liminary work by Spurzem & Einsel (1998). In this paper, 
thus we present the study of dynamical evolution of the ro- 
tating star cluster with the initial mass function (IMF here- 
after). We first investigate the evolution of two-component 
models as a simplest extension of the single mass models. 
Then we extend our study to the multi-mass models repre- 
sented by ton different mass species. The difference of our 
models as compared to the preliminary work by Spurzem & 
Einsel (1998) is that first our improved code includes more 
accurate numerical integration and discretisation procedures 
as described in Paper II, and here we do an extensive pa- 
rameter study, which did not exist before. 

The outline of the present paper is as follows; In section 
2, the FP equations for the multi-mass system and the initial 
models are presented. We concentrate on the evolution up 



to the core collapse of two component models in section 3 
and multi-mass models in section 4. Wc further discuss the 
evolution beyond the core collapse for both two-component 
and multi-component models section 5. We summarize our 
main results in section 6. 



2 THE MODELS 

2.1 Fokker-Planck Equations 

We have constructed a computational scheme to study the 
dynamical evolution of rotating stellar systems with mass 
spectrums in detail with higher accuracy both for pre- and 
post-collapse. The framework of the method is essentially 
an extension of the method used in Paper II. The multi- 
mass FP equation under a fixed potential can be written in 
a flux-conserving form as follows: 



P(E J)^--^ 



dFj 



(1) 



where P{E,Jz) is the phase space volume accessible for 
stars with E and Jz, and and Fj^ ^ are the dis- 

tribution function and the particle flux in energy {E) and 
a-component of angular momentum (Jz) of the i-th compo- 
nent, respectively. The expression for phase space volume is 
given in Paper I. The particle fluxes Fe^ and Fj^ . can be 
expressed as follows: 



dfi dfi 

-Fsf = DEE,i-^ + DEJ^,i^r-r- + DE,ifi; 

O-ty OJz 



-Fj^ . = Dj,j,,i^ + Dj,E,i^ Dj^ Ji. 
dJz oE 



(2) 



where Dee.x, DEj^,i, etc are the diffusion coefficients and 
are given in the Appendix of Paper I for single component 
models. The extention to the multi-mass models can be sim- 
ply done by applying the distribution function of the i-th 
component in the expression of these coefficients. 

It is necessary to add energy source to explore the evo- 
lution after core-collapse. Primordial binaries and massive 
stars are important energy sources, since they can delay 
significantly the core-collapse time and change the early 
evolution of stellar systems (Gao et al. 1991, Giersz & 
Spurzem 2000). However, the specification of primordial bi- 
naries could introduce many more model parameters. The 
treatment of such binaries in FP code is not so trivial ei- 
ther. Since we are mainly concerned with the general effects 
of rotation on the dynamical evolution of stax clusters, we 
limit ourselves to simplest cases. As a way to follow the dy- 
namical evolution beyond the core collapse, wc have only 
considered the heating effect of binaries formed by three- 
body processes. Heating formula for three-body binaries of 
Lee, Fahlman & Richer (1991) is used in our models (see 
also Takahashi 1997 for the discussion of this formula). In 
their formulation, the total heating rate per unit volume is 
given by 



2 



(3) 
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where nj, rrii, and ai are the number density, stellar mass, 
and one-dimensional velocity dispersion, respectively, of the 
i-th component, and ao is the mass density-weighted central 
one-dimensional velocity dispersion. The standard value of 
Cb ~ 90 is used for present study. Heating rate for each mass 
component is obtained by distributing the total heating rate 
such that 

Etot.i = -^—Etot (4) 

Ptot 

where pi and ptot represent mass density of each component 
and total density. Gicrsz & Spurzem (1997) have demon- 
strated that such a simple approximations work very well 
for modelling the post-collapse evolution of globular clus- 
ters. 



2.2 Multi-Mass Initial Models 

Wo employ the rotating King models as initial models fol- 
lowing Lupton & Gunn (1987). We assume that there is 
no initial mass segregation among different mass species. 
Our present initial models are, therefore simple extension 
of the initial models used in Papers I and II. These models 
arc characterized by two parameters: dimcnsionlcss central 
potential Wo and rotational parameter uiq. We have stud- 
ied the evolution of clusters with Wq = 6 and Wo = 3, 
respectively. Model clusters are assumed to rotate in circu- 
lar orbits around mother galaxy, so that the mean density 
within the tidal radius (rt) remains a constant throughout 
the evolution. We examined the evolution of clusters with 
various amount of the initial rotations (see Tables 1 and 
2 for the list of the initial models of two component, and 
multi-component, respectively). 

First, we considered the models with two-component 
mass species. Two-component models are expected to pro- 
vide us with understanding the essential features of the evo- 
lution of the multi-component clusters (for example, the 
transfer of angular momentum between two mass compo- 
nents). We have examined the wide range of individual stel- 
lar mass ratio fx := m2/rn\ and total mass ratio M1/M2 = 
(1 - q)/q, where q := M2/(Mi -|- M2). The quantities n 
and q can be compared with the work of Khalisi (2002) and 
Fregeau et al. (2002). 

Next, we considered clusters which have continuous 
mass spectra. We choose simple power-law mass function 
for convenience. The number of stars in a mass interval 
(m, m + dm) is given by 

dN{m) = Cm"dm, mmin < m < mmax (5) 

where mmin and mmax denote the minimum and the maxi- 
mum masses, respectively. The dynamic range of mass spec- 
trum, expressed as mmax/mmin is 10 for present study. 

We have studied the evolution of clusters with three 
different shapes of mass function: a— —1.20, —2.35 and 
—3.50, where a = —2.35 represents Salpeter-type mass func- 
tion. We list the parameters of the models with continuous 
meiss spectrum in Tables 1 and 2. For models of post core- 
collapse, we consider only a model with a mass function with 
a = -2.35. 

The number of mass groups used in present study is 
ten. The mass of each group is assigned such that 



Table 1. Initial Models of clusters with Wo = 6. 



Model 1712 /mi M1/M2 


ao 




[%] 


Phase 






0.0 


0.00 


Pre 






0.2 


3.40 


Pre 






0.3 


7.06 


Pre 






0.4 


11.35 


Pre 






0.6 


20.17 


Pre 


M2A 2 5 


— 


0.8 


27.96 


Pre 






1.0 


34.23 


Pre 






1.2 


39.13 


Pre 






1.4 


42.93 


Pre 






1.6 


45.88 


Pre 






0.0 


0.00 


Pre 


M2B 2 10 





0.3 


7.06 


Pre 






0.6 


20.17 


Pre 






0.0 


0.00 


Post 


M2C 5 in 


_ 


n.3 


7.06 


Post 






0.() 


21). 17 


I'ost 






0.0 


0.00 


Pre 


M2D 5 100 





0.3 


7.06 


Pre 






0.6 


20.17 


Pre 






0.0 


0.00 


Pre 


M2E 10 10 





0.3 


7.06 


Pre 






0.6 


20.17 


Pre 






0.0 


0.00 


Pre 


M2F 10 100 





0.3 


7.06 


Pre 






0.6 


20.17 


Pre 






0.0 


0.00 


Pre 


MCA 


-1.20 


0.3 


7.06 


Pre 






0.6 


20.17 


Pre 






0.0 


0.00 


Post 






U.o 


* .UD 


Post 






0.6 


20.17 


Post 






0.0 


0.00 


Pre 


MCC 


-3.50 


0.3 


7.06 


Pre 






0.6 


20.17 


Pre 


Tr'. total rotational energy 










Ti^: total kinetic energy 










rui = mmin 

\ mmin / 


, (i = 1, 


2,.. 


.10). 


(6) 


The total mass of ith mass group is 








»",■ I 1 
2 










Mi= / N(m)mdm. 








(7) 



m . 1 

If the initial density distribution and the initial mass 
function including mmin and m^nax arc fixed, the only free 
parameter for the model is the number of stars (or total 
mass if the masses of stars are given in physical units). The 
general behavior of the core collapse docs not depend on 
the number of stars. However, the core stops to collapse and 
begins to expand when the central density exceeds a cer- 
tain value which depends on the total number of stars in 
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Table 2. Initial Models of clusters with Wo = 3. 







Ml /Mo 




UJQ 


T IT, fO"! 














[%] 












0.0 


0.00 


Pre 










0.4 


1.89 


Pre 


Aifo A 
MzA 


2 


5 




0.8 


6.96 


Pre 










1.2 


13.84 


Pre 










1.5 


19.37 


Pre 










1.6 


21.20 


Pre 










0.0 


0.00 


Pre 


M2B 


2 


10 





0.8 


6.96 


Pre 










1.5 


19.37 


Pre 










0.0 


0.00 


Post 


M2C 


5 


10 





0.8 


6.96 


Post 










1.5 


19.37 


Post 










0.0 


0.00 


Pre 


M2D 


5 


100 





0.8 


6.96 


Pre 










1.5 


18.37 


Pre 










U.U 


U.OU 


Pre 


M2E 


10 


10 





0.8 


6.96 


Pre 










1.5 


19.37 


Pre 










0.0 


0.00 


Pre 


M2F 


10 


100 




0.8 


6.96 


Pre 










1.5 


l<).:-!7 


I'rc 










0.0 


0.00 


Pre 


MCA 






— 1.20 


0.8 


6.96 


Pre 










1.5 


19.37 


Pre 










0.0 


0.00 


Post 


MCB 






-2.35 


0.8 


6.96 


Post 










1.5 


19.37 


Post 










0.0 


0.00 


Pre 


MCC 






-3.50 


0.8 


6.96 


Pre 










1.5 


19.37 


Pre 



core bounce due to the formation of binaries. Aside from 
the actual values of the central density, the general behavior 
of the post-collapse evolution does not sensitively depend 
on initial N. Note that there will be rapid core oscillation, 
which is beyond the scope of the present study, if we use 
realistic particle numbers. 



2.3 Tidal Boundary 

All models studied in the present work are tidally limited. 
Realistic treatment of tidal field is possible only for direct 
N-body calculations. The simplest assumption for the tidal 
boundary is to assume a spherical shape whose radius is 
determined by the requirement of constant mean density 
within this radius. The actual shape is not exactly spherical 
and the tidal field varies with time as the cluster moves 
along non-circular orbit in non-spherical Galactic potential. 
Nevertheless, we assumed a constant mean density within 
the tidal radius for simplicity. 

We have applied the "energy condition" for the removal 
of stars (see Takahashi, Lee & Inagaki 1997 for details of 
tidal conditions). Takahashi & Portegies Zwart (1998) found 
that the "apocenter condition" reproduces the result of the 
N-body simulation better for small N models, but the differ- 
ence between different conditions diminishes as A'^ becomes 
larger. Since our distribution functions do not depend on 
total angular momentum, we could not apply the apocen- 
ter condition. However, we emphasize that our choice of 
N = 10, 000 is simply to reduce the computational time. 
Although the parameters for the central parts (central den- 
sity, core radius and central velocity dispersion) are sensitive 
to the N , the global evolution (such as M(t), rh{t)) do not 
depend on AT as long as we use laxge AT approximation for 
the tidal conditions. We could consider our results to be ap- 
propriate for systems with very large N , except for the core 
properties. 



the cluster. If N is too small, we cannot apply the Fokker- 
Planck method to follow the evolution of the star cluster. 
The development of gravothermal oscillation occurs for N 
much larger than 10,000 because the post-collapse cluster 
becomes gravothcrmally unstable (Goodman 1987). The use 
of larger N also takes more computational time because the 
post-collapse expansion begins at much higher central den- 
sity. As a compromise, we have used N = 10, 000 for models 
in the present study. 

The number of stars used here is quite small compared 
to that of globular clusters. According to the recent compar- 
isons between direct A''-body models and the Fokker-Planck 
models (Takahashi & Portegies Zwart 1998, Portegies Zwart 
& Takahashi 1999) the evolution of the total mass until total 
dissolution depends on the particle number since the mass 
loss process involves the crossing time rather than the relax- 
ation time. Nevertheless we restrict ourselves here to post- 
collapse models using relatively small N , because this speeds 
up the computational time, and we are mainly interested in 
the interplay between rotation, relaxation and tidal mass 
loss. N = 10, 000 is large enough to study the effect of the 



3 PRE-COLLAPSE EVOLUTION OF 
TWO-COMPONENT MODELS 

Under the context of our models of two-body relaxation, ro- 
tation, and binary heating the clusters always go through the 
corc-coUapsc and l)cgins the post-collapse expansion. Some 
aspects of the evolution differ between pre- and post-collapse 
phases. Thus we first discuss the pre-coUapse evolution of 
two-component models first. In the study of the prc-collapse 
phase only, we have ignored the heating by binaries. 

For clusters with multiple mass, the energy exchange 
between different mass species causes mass segregation 
which takes place in a time scale shorter than the single- 
mass core-collapse. Once the core is dominated by the mas- 
sive component whose velocity dispersion is lower than aver- 
age, the core collapse proceeds more rapidly than the single 
component clusters. The rotation is also known to accelerate 
the core-collapse process since it provides additional mecha- 
nism of driving the catastrophic collapse. In two-component 
clusters, these two processes compete with each other. We 
now investigate how the rotation affects the course of core- 
collapse of two-component models. 



© 0000 RAS, MNRAS 000, 000-000 



Dynamical Evolution of Rotating Stellar Systems: III. The Effect of Mass Spectrum 5 




Tidally limited, M,/Mj=5 
I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 




0.2 0.4 0.6 O.B 1 1.2 1.4 1.6 



Figure 1. Time evolution of the central density for the cluster 
model with mass function M2A 

3.1 Central density and core-collapse time 

Wo have shown the evolution of central density for the model 
M2A with different initial central potential, Wo = 3 and 6 
in Fig. 1. he time is expressed in units of the initial half- 
mass relajcation imie{Trhfi), and the central density in units 
of Mo/r^Q, where Mo and rc,o are the initial mass of the 
cluster and the initial King's core radius, respectively. An 
expression for the initial half-mass relaxation time is taken 
from Spitzer & Hart (1971), 

_(U38iV^^ 

^-".o - Gi/2^i/2inA 

where No,rh,o,G,fh and In A denote the initial total number 
of stars the initial half-mass radius, gravitational constant, 
the moan mass of the particles and the Coulomb logarithm, 
respectively. Since the rotating clusters are not spherically 
symmetric, some cautions should be taken in defining the 
half-mass radius. It is defined as the effective radius of con- 
centric spheroid within which the enclosed mass becomes 
half of the total mass (see eq. [22] of Einsel & Spurzem 
1999). 

It is shown clearly that the clusters with higher initial 
rotation reach core-collapse earlier than the cluster without 
initial rotation. In order to quantify the amount of accelera- 
tion of the evolution, we define the "degree of acceleration" 
as follows: 

^ tccjui = 0) - tccjuj) 

dace — J- f \ ' y'^l 

where tec is the time to reach the core-collapse. Obviously, 
dace depends on u> for a given Wo and a mass function. We 
have listed dace in the last column of Tables 3 and 4, and 
have shown them in Fig. 2(b). 

Fig. 2(a) shows the run of the core-collapse times as 
a function of the initial rotation (solid line) and the ratio 



Figure 2. (a) Run of core coUpase time measured in initial half 
mass rclexation time as a function of the initial rotation rotational 
parameter (ioo) for a scries of models with a mass function M2A 
and Wo = 6 (solid line). A ratio between the initial rotational 
energy and the total kinetic energy for the same models (dashed 
line), (b) The degree of the acceleration of core-collapse defined 
by eq. 9 for model with the central potential Wo = 6 (solid line) 
and 3 (dashed line), respectively. 



between initial total rotational energy due to the rotation 
and initial total kinetic energy (dashed line) for model with 
mass function of M2A and the initial central potential Wo = 
6. The ratio between rotational energy and kinetic energy 
of the model with initial rotation wo = 1-6 is ~ 46%, just 
below the dynamically unstable criterion of Tr/Tk = 0.5, 
where Tr and Tk denote the total rotational energy and total 
kinetic energy. The coro-coUapsc times, which arc measured 
in units of the initial half-mass relaxation time decrease as 
the initial degree of rotation increases. The cluster with the 
initial rotation uo = 0.6 reaches the core collapse faster by a 
factor of ~ 1.3 than the cluster without the initial rotation. 
The cluster with the initial rotation ujq = 1.6, a model with 
the largest initial degree of rotation among current models, 
has a core-coUase time of just ~ 1/3 of that of the model 
without the initial rotation. 

We have shown the time evolution of central density for 
models with different mass function except for mass function 
of M2A in Fig. 3. Solid lines represent the models without 
initial rotation. Dotted lines and dashed lines represent the 
models with the initial degree of rotation ujq — 0.3 and 0.6, 
respectively. We have chosen the values for the initial degree 
of rotation as the same values used in Paper II. The mass 
ratios between high and low mass components(m2/mi) and 
the ratios of total masses of these components (M1/M2) are 
written to distinguish the different mass functions. 

In general the core collapse occurs earlier for more 
rapidly rotating models. The degrees of acceleration also 
depends sensitively on the assumed mass function. For ex- 
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Figure 3. Time evolution of the central density for the cluster 
model with two-component mass functions and Wo = 6 except 
for a mass function of M2A. 

ample, dace of the model with the initial rotation uio = 0.6 
and mass function of M2B is much larger than that of the 
model with mass function M2E, which has the same initial 
degree of rotation. The main difference between these two 
models is the ratio of the individual mass: m2/mi = 2 for 
M2B and 10 for M2E. The larger difference in individual 
mass means larger amount of energy exchange and thus the 
significant amount of acceleration of core collapse due to 
this effect alone. There appears not enough room for the 
acceleration due to the initial rotation for such cases. For 
a given 1712 /mi, the degree of acceleration also depends on 
M1/M2. If the total mass in high mass component is rela- 
tively large (i.e., smaller values for M\/M2), the acceleration 
due to rotation is small. Such trend can be seen from the 
comparison between M2C and M2D or M2E and M2F. For 
models with relatively small M1/M2, the acceleration due 
to mass segregation is already significant and the role of 
rotation is relatively less important. 

Time evolution of total mass of the clusters with mass 
function M2A is displayed in Fig. 4. The total mass of each 
mass component is normalized by their initial total mass 
Mifi. The dashed lines represent the total mass of lighter 
component and the solid lines for massive component, re- 
spectively. The total mass of lighter component decreases 
more rapidly than the that of massive stars, irrespective of 
the initial degree of rotations. The cluster with a initially 
higher rotation loses mass more efficiently than the clusters 
with lower initial rotation. The higher moss-loss rate of clus- 
ter with the higher initial rotation was known in previous 
works for single mass models (Papers I and II). The ratio 
of total masses at a time of core-collapse between two mass 
components varies with different initial degree of rotations. 
For the model with initial rotation of wo = 1-6, the total 
mass of the low mass component at the time of core-collapse 



Figure 4. Evolution the total masses retained in a cluster for 
model with Wa = 6 and mass function M2A. The total mass of 
the individual component is normalized by its initial total mass 
Mifi. Solid lines represent the evolution for the massive stars and 
dashed lines for the less massive stars, respectively. 

is ^ 10 % of its initial value. On the other hand, ~ 30 % 
of the initial mass still remains for massive component. The 
ratio of total mass at a time of core-collapse between the 
high and low mass components for model without the initial 
rotation is ~ 1.1, which is much lower than the value ob- 
tained for cluster model with the initial rotation wo = 1.6. 
It implies that the mass function of rotating cluster which 
has the same initial mass function will follow different evo- 
lution according to their initial rotation. The evolution of 
mass function, especially the power law index q, for models 
with initially power-law mass function is explained in detail 
in Sect. 5.4. 



3.2 Central velocity dispersion and central 
anguleir speed 

Einsel & Spurzem (1999) argued in their single mass model 

that the acceleration of corc-coUapse time for models with 
the initial rotation compared to the model without the initial 
rotation is accompanied by a rapid increase of the central 
velocity dispersion in rotating models. We have checked if 
this is also true for two component models. Fig. 5 shows the 
evolution of mass density weighted central one-dimensional 
velocity dispersion (ctc) for models with the central potential 
W o = 6 and mass function model of M2A. The initial rota- 
tional parameters (cJo) are written in front of the starting 
points of evolution. 

The evolution of Gc measured in the total central den- 
sity (pc) shows a rather complex behavior. For models with 
large initial rotation, ctc increases in the very early phase of 
evolution and then decrease until it reaches broad minimum. 
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Figure 5. Evolution of the mass-weighted central ID velocity 
dispersion with respect to the central density for cluster models 
with Wtt = 6 and M2A for a mass function. Tlic values of initial 
rotational parameter (ujo) are written ahead of the starting points 
of the evolution. 



Figure 6. Evolution of the central angular speed of the individual 
mass species with respect to the total central density for cluster 
models with H-o = 6 and M2A for a mass function. The central 
rotational speed of the massive stars (solid lines) increases very 
fast during the early evolutionary stage, while that of the lighter 
stars (dashed lines) shows an moderate increase. 



Eventually, ctc increases with pc monotonically toward the 
core collapse. 

The initial rise of Uc is caused by the initial collapse and 
subsequent heating by contraction. The decrease of (Jc after 
the local maximum is due to the loss of energy of high mass 
stars to low mass stars. Since there are more low mass stars 
in the early phase in the central parts, the mass-weighted 
velocity dispersion decreases. The subsequent increase of o"c 
after the broad minimum occurs when the central part is 
dominated by high mass stars. The stars in the central parts 
lose energy to stars in the outer parts and the velocity dis- 
persion increases because of the negative specific heat. 

The evolutionary path reaches self-similar collapse 
phase when the central density increases by 2 ~ 3 orders 
of magnitude. The evolution of self-similar collapse can be 
characterized by a parameter given by 



7 = 



dlncr^ 
din Pc 



(10) 



For single mass system, Einsel & Spurzem (1999) obtained 
an average value of 7 = 0.109, independent of the initial de- 
gree of rotation. It is expected that 7 for present two compo- 
nent models is very similar to that of the single component 
system, since the central part of cluster is dominated by high 
mass stars. We derive an average value of 7 = 0.105, which 
is very similar to what was obtained for single mass model. 

Another feature of Fig. 5 is the initial evolution of of 
(Tc depends on the amount of initial rotation. The minima 
occur at higher central density for models with larger initial 
rotation. The enhanced two-body relaxation process due to 
rotation will eventually cause more rapid increase the central 



density. This results in the acceleration of the core-collapse 
for rapidly rotating systems. 



Fig. 6 displays the run of central angular speed (f2c) as 
a function of the total central density for the model with 
the mass function of M2A. The solid lines and the dashed 
lines represent fie of high mass and low mass components, 
respectively. It is shown clearly that f2c of the higher mass 
component increases very quickly only during early phase of 
evolution. For the models rotating slowly (wo = 0.2 and 0.4) 
fie of lower mass component decreases by a small amount. 
For rapidly rotating clusters, however, f7c increases with in- 
crease of central density. The rapid increase of f2c of high 
mass component in the early phase is a caused by rapid col- 
lapse of the high mass stars: as they lose energy, they rotate 
faster. 

The mass weighted angular speeds in the core as a func- 
tion of Pc for the models in Fig. 6 are shown in Fig. 7. During 
early evolution, f2c increases rapidly for models with large 
rotation, then it grows slowly with increase of central den- 
sity, reaching self-similar phase in evolutionary path. Einsel 
& Spurzem (1999) introduced the parameter <5 which char- 
acterizes the self-similar evolution of central angular speed 
as follows, 



5 = 



din f2c 



(11) 



din pc 

The calculated values of 8 for the present two- 
component models are listed in Table 3. On average, <5 sa 
0.110 with weak dependence on the initial rotation for model 
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Table 3. Physical properties of cluster models with Wo = 6. 



ivioaei 


CJO 


tCc/Trhfl 




tdis/Trh,0 


ctacc 


7 


R 





0.0 


5.74 


0.81 




0.00 


0.108 






0.2 


5.57 


0.78 


- 


0.03 


0.114 


0.112 




0.3 


5.36 


0.72 


— 


0.07 


0.106 


0.114 




0.4 


5.06 


0.66 


— 


0.14 


0.106 


0.113 


M2A 


0.6 


4.33 


0.53 




0.33 


0.105 


0.111 




0.8 


3.54 


0.41 





0.62 


0.105 


0.108 




1.0 


2.83 


0.31 


_ 


1.03 


0.105 


0.108 




1.2 


2.30 


0.24 




1.50 


0.105 


0.111 




1.4 


1.92 


0.18 




1.99 


0.105 


0.109 




1.6 


1.76 


0.13 




2.26 


0.106 


0.107 




0.0 


6.69 


0.78 




0.00 


0.106 




M2B 


0.3 


6.13 


0.70 




0.09 


0.105 


0.108 




0.6 


4.76 


0.52 




0.40 


0.104 


0.111 




0.0 


1.28 


0.96 


10.68 


0.00 


0.106 




M2C 


0.3 


1.27 


0.93 


7.74 


0.01 


0.105 


0.090 




0.6 


1.22 


0.82 


4.36 


0.05 


0.104 


0.092 





0.0 


4.29 


0.88 




0.00 


0.101 




M2D 


0.3 


3.96 


0.82 




0.08 


0.099 


0.092 




0.6 


3.13 


0.67 




0.37 


0.102 


0.094 




0.0 


0.48 


0.98 




0.00 


0.106 




M2E 


0.3 


0.47 


0.97 




<0.01 


0.105 


0.090 




0.6 


0.47 


0.90 




0.01 


0.104 


0.084 




0.0 


0.84 


0.98 




0.00 


0.099 




M2F 


0.3 


0.80 


0.96 




0.05 


0.099 


0.099 




0.6 


0.77 


0.90 




0.08 


0.101 


0.100 




0.0 


3.49 


0.83 




0.00 


0.106 




MCA 


0.3 


3.39 


0.74 




0.03 


0.104 


0.095 




0.6 


2.83 


0.54 




0.23 


0.103 


0.099 




0.0 


1.98 


0.92 


9.60 


0.00 


0.104 




MCB 


0.3 


1.92 


0.86 


7.11 


0.03 


0.103 


0.105 




0.6 


1.72 


0.70 


4.16 


0.15 


0.103 


0.113 




0.0 


1.73 


0.94 




0.00 


0.098 




MCC 


0.3 


1.68 


0.90 




0.03 


0.096 


0.105 




0.6 


1.55 


0.77 




0.12 


0.096 


0.103 



tec /'''rhfi'- core collapse time measured in the initial half-mass relaxation time 
Mqc /Mq: mass of the cluster in unit of the initial total mass (Mq) at t = tec 
t^is/r^hfi- time in unit of r^^fi when the cluster dissolved completely 
7, <5: see the text for the definitions. 



M2A. This value of 5 is about twice of those of single mass 
models. Einsel & Spurzem (1999) obtained 6 = 0.06 ~ 0.08 
depending on the initial rotation for their model cluster with 
central potential Wo = 6. 

The relationship between the central density, ctc, and 
Q,c for models with mass function M2E (m2/mi = 
10, M1/M2 = 10) are shown in Fig. 8. For Uc versus pc, wc 
have plotted the velocity dispersion of individual mass com- 
ponent for model without initial rotation and models with 
Wo = 0.3 and 0.6, while only the rotating models are shown 
for Oc versus pc plot. The starry mark in Fig. 8(a) represents 
the starting point of evolution for both mass components. 
The central velocity dispersion of high mass component de- 
creases dramatically, while that of low mass component in- 
creases slightly at early evolutionary stage due to the energy 



cquipartition between two mass components. The decrease 
of the central velocity dispersion for model with higher ini- 
tial rotation is larger than that of the lower initial rotation. 
Inspection of Fig. 8(a) reveals that the slope in log- log plot 
for massive component becomes 7 « 0.105 regardless of de- 
gree of rotation. 

3.3 Rotational Velocity and Angular Momentum 

Two-body relaxation causes the outward transfer of the an- 
gular momentum and escape of stars at the outer parts, 
leading to the decrease of rotational energy. Kim ct al 
(2002) showed that the majcimum value of rotational veloc- 
ity (Kot.max) in equator decreases monotonically for their 
single mass models. In Fig. 9, we displayed the evolution of 
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Table 4. Initial Models of clusters with Wo = 3. 



Model 






Mcc/Mq 




dace 


7 


<5 




0.0 


5.78 


0.42 


- 


0.00 


0.107 


- 




0.4 


5.60 


0.39 


— 


0.03 


0.106 


0.110 


MzA 


0.8 


5.05 


0.33 




0.15 


0.106 


0.112 




1.2 


4.27 


0.25 




0.35 


0.106 


0.113 




1.5 


3.61 


0.18 




0.60 


0.106 


0.109 




1.6 


3.37 


0.15 


- 


0.71 


0.106 


0.111 




0.0 


6.21 


0.41 




0.00 


0.105 




M2B 


0.8 


5.37 


0.33 





0.16 


0.104 


0.107 




1.5 


3.79 


0.19 


- 


0.64 


0.104 


0.113 




0.0 


1.98 


0.76 


4.17 


0.00 


0.105 




M2C 


0.8 


1.88 


0.66 


3.43 


0.05 


0.104 


0.094 




1.5 


1.59 


0.39 


2.01 


0.25 


0.105 


0.090 




U.U 


A OO 
4.0Z 


U.o7 




U.OU 


0.104 




M2D 


0.8 


4.82 


0.49 





0.18 


0.103 


0.095 




1.5 


2.80 


0.34 


- 


0.72 


0.104 


0.099 




0.0 


0.86 


0.87 




0.00 


0.105 




M2E 


0.8 


0.84 


0.80 




0.02 


0.104 


0.084 




1.5 


0.76 


0.56 


- 


0.13 


0.104 


0.085 




0.0 


1.54 


0.86 




0.00 


0.102 




M2F 


0.8 


1.49 


0.79 





0.03 


0.101 


0.087 




1.5 


1.37 


0.60 


- 


0.13 


0.101 


0.093 




0.0 


4.01 


0.38 




0.00 


0.104 




MCA 


0.8 


3.51 


0.30 




0.14 


0.105 


0.106 




1.5 


2.51 


0.15 




0.60 


0.105 


0.103 




0.0 


2.65 


0.56 


4.18 


0.00 


0.102 




MCB 


0.8 


2.42 


0.45 


3.48 


0.09 


0.101 


0.142 




1.5 


1.83 


0.23 


2.31 


0.45 


0.103 


0.139 




0.0 


2.48 


0.70 




0.00 


0.096 




MCC 


0.8 


2.29 


0.58 




0.08 


0.096 


0.118 




1.5 


1.82 


0.34 




0.36 


0.096 


0.135 



maximum rotational speed at equator, where the effect of 
the rotation is largest, for model with the central potential 
Wo = 6 and mass function M2A. Only four different models 
for initial degree of rotation is plotted to avoid complexity. 
Different line styles represent the different initial rotations. 
The thick lines represent the run of maximum rotational 
velocity at equator for the high mass component, and the 
thin lines for low mass component. The maximum rotational 
speed of low mass component decreases with the central den- 
sity. The central density increases with time monotonically 
for core-collapse models as shown in Sect. 3.1. The rotational 
speed for the high mass component, especially for model 
with the initial degree of rotation ujo = 0.3 increases during 
early evolutionary stage, and it remains around 0.2 after- 
ward. However, the maximum rotational speed for rapidly 
rotating model (i.e., uJo = 1.6) decreases with central den- 
sity monotonically for both components. We have shown 
the run of the maximum rotational speed at equator for 
models with the central potential Wo = 6 and too = 0.6 
in Fig. 10. To investigate the difference in rate of angular 
momentum transfer between two mass components, we in- 
clude 6 different mass function models. The individual mass 



ratio (mg/mi) and the total mass ratio (M1/M2) of initial 
model are marked in the figure. While the maximum rotar 
tion speeds of the high mass component increases slightly 
for model with lowest individual mass ratio (m2/mi — 2), 
it increase more steeply for other models. The maximum ro- 
tational velocity of high mass component for models with 
m^/mi = 5 and 10 remains a nearly constant afterward. 



Kim et al. (2002) found that the specific angular mo- 
mentum (Jz,max) decreases with time for models with cen- 
tral King's potential Wo = 6. Fig. 11 shows the evolution 
of maximum specific angular momentum with the increase 
of central density for the same models to those of Fig. 10. 
The line style and the thickness of each line represent the 
same mass function model shown in Fig. 10. The maximum 
angular momentum decreases monotonically irrespective of 
the mass function and the mass component, although the 
rate of decrease depends on the individual mass ratio and 
the mass components. 
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Figure 7. Evolution of the mass-weiglited central angular speed 
with respect to the central density for cluster models with Wo = 6 
and M2A for a mass function. 



W„=6, Tidally limited, 1112/111,= 10, M,/M2=10 




0.1 1 10 100 1000 10< 10= 





(b) 














— ... - 0.3 : 












solid; m = mg | 
dashed; m^rrij | 



0.1 1 10 100 1000 lO" 

Pc 



Figure 8. Evolution of the central velocity dispersion (a) and the 
central angular speed (b) for the models with Wo = 6 and a mass 
function M2E (m2/mi = 10, M1/M2 = 10), respectively. The 
start of the evolution is marked with starry mark in figure (a). 
We have included models with the initial rotation of uo = 0.0, 0.3 
and 0.6. However, only the rotating models are considered in Qc 
versus pc plot. The steep decreasing of the central velocity disper- 
sion and the steep rising for central angular speed of the massive 
component are shown clearly. The loss of central velocity disper- 
sion is compensated by the quick increasing of angular speed. 



Figure 9. Evolution of the maximum rotational velocity at the 
equator according to the central density for models with Wo = 6 
and mass function of M2A. We have shown only for four different 
initial rotations out of 9 models. The thin lines and the thick 
lines represent the evolution of the lighter component and that of 
massive component, respectively. 



0.4 



e 0.3 



0.2 



0.1 1 10 10^ 10= lO* 10= 

Pc 

Figure 10. Same as the Fig. 9 but for all two-component models 
with uiQ = 0.6 of present study. The Vrot,max of the massive stars 
for the models with higher individual mass ratio (m2/mi) even 
increases very fast during the early evolutionary stage due to the 
transfer of the angular momentum from the lighter stars. 
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(W(,,M„) = (6,0.6), Tidally limited, 2-comp. models Tidally limited, 10-comp. models 




Figure 11. Evolution of the maximum specific angular momen- 
tum at the equator for the cluster model with two-component 
mass functions. The meaning of the line styles are the same to 
those of Fig. 3-10. 

4 CLUSTERS WITH CONTINUOUS MASS 
SPECTRUMS 

4.1 Evolution of central properties 

Fig. 12(a) and (b) shows the time evolution of the total 
central density for cluster models with mass functions MCA, 
MCB and MCC. For each mass function three models with 
the different initial rotation are plotted together. Since we 
select the continuous mass function of the power-law type 
and the index of power-law varies with time due to difference 
in mass evaporation rates among different mass groups, we 
distinguish the initial power-law index (qq) with the power- 
law index (a) at time t. The values of the initial power-law 
index are shown in figure and the models with the same 
mass function are connected with lines for clarity. 

Like single component or two-component models, the 
models with the higher initial rotation evolve more rapidly 
than the those with lower initial rotation. Among the models 
with the same amount of the initial rotation the cluster with 
the flatter initial mass function (qq — —1.20) has a longer 
collapse time, when the time is measured in the initial half- 
mass relaxation time. For models without initial rotation it 
is known that the core-collapse time for cluster with steeper 
meiss function is shorter than that for the cluster with shal- 
lower mass function (e.g., Lee, Fahlman & Richer 1991; Lee 
& Goodman 1995; Takahashi 1997). 

We have shown the evolution of the central one dimen- 
sional velocity dispersions and the central angular speeds of 
individual mass components for the with {W»,ljo) = (6, 0.6) 
and the initial power-law slope qq = —2.35 in Fig. 13. The 
evolution starts at lower left corner as marked by a starry 
point in each figure. The individual mass of each mass group 
decreases upward in Uc versus pc, and vice versa for f2c ver- 



Figure 12. Time evolution of the central density for the cluster 
model with the central concentration Wo = 6 (a) and 3 (b) and 
contiimous mass function (10 mass components). Three different 
initial rotation are considered. The cluster models with the same 
mass function and different initial rotation is connected with a 
line. The acceleration of the core-collapse due to the initial rota- 
tion is shown clearly for all models. 



sus pc as shown in Fig. 13. The parameters 7 and 5 which 
characterize the core-collapse phase as defined in Eqs. (10) 
and (11) are 7 = 0.103 and 6 = 0.113. 

The effects of the initial mass function on the central 
velocity dispersion and the central angular speed are shown 
in Figs. 14 15, respectively. We keep the scales of the hori- 
zontal axis and the vertical ax;is to be the same for all panels. 
The direction for the increase of mass of each mass group is 
indicated by arrows in the figures. For most massive com- 
ponent as ao decreases Uc decreases more rapidly, driving 
more rapid evolution of the cluster as indicated by earlier 
core-collapse. The development of the mass segregation is es- 
tablished more quickly for the model with the steeper mass 
function. For Oc versus pc (Fig. 15) we have shown only the 
early stage of evolution since the rotation affects mainly dur- 
ing the initial period. It is evident that the development of 
self-similax core-collapse phase occurs earlier for the model 
with the steeper meiss function than that with fiatter mass 
function. The central angular speed at a given central den- 
sity is larger for the model with the steeper mass function. 
The rapid increase of the central angular speed for the high- 
est mass group and the slow decrease for the lowest mass 
group for the model with ao = —3.50 is a consequence of 
the rapid development of mass segregation. The steep de- 
crease of central velocity dispersion drives the increase of 
the central angular speed because of negative specific mo- 
ment of inertia. 
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(W„,cj„) = (6,0.6) Ticlally limited, a„- a.'iS 




0.1 1 10 10= 10' lO* 



Figure 13. Evolution of the ID central velocity dispersion (a) 
and the central angular speed (b) of the individual mass compo- 
nents for a cluster model with (l-Vi), wo) = (6, 0.6) and a power- law 
mass fimction with index ao = —2.35. The staray marks denote 
the start of evolution. Not only the velocity dispersion, but the 
central angular speed shows a self-similar behaviour. 



4.2 Rotational properties 

Fig. 17 displays the evolution of the maximum rotational 
velocities of individual mass components at the equator for 
the model with {Wo,uJo) = (6,0.3). We included the mod- 
els with mass functions MCA, MCB and MCC together. 
The total mass of each mass group increases along the ax- 
row shown in Fig. 17(a). The evolution of rotational speed 
for the cluster with the continuous mass function is simi- 
lar to that of the two component models. For the model 
with ao = —3.50 the maximum rotational speed of the mass 
component of ms remains nearly a constant. Since the ra- 
dius where the rotational velocity becomes maximum at the 
equator decreases with the time, the mass group with ms 
rotates faster with time.The constant behaviour of the max- 
imum rotational speed occurs for the mass group of my for 
model with ao = —2.35 and mass group mg for model with 
Qo — —1.20, respectively. Fig. 16(a) and (b) show the evo- 
lution of the maximum rotational speed and the evolution 
of the maximum angular momentum at equator for model 
with (Wo, Wo) = (6, 0.6) and mass function with ao = —2.35. 
The maximum rotational speed of the lower mass groups 
decreases more than the cluster rotating slowly (see Fig. 
17(b) for comparison). While the high mass component ro- 
tates faster, the angular momentum decreases continuously 
through whole evolutionary phase. Since the outward trans- 
fer of angular momentum for each mass species is closely 
related with the two-body relaxation process, the evolution 
of the maximum of the angular momentum reaches a plateau 
when it is measured in the central density. 



(W„,a;o) = (6,0.3), Tidally limited, 10-comp. models 
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Figure 14. Evolution of the ID central velocity dispersion for 
cluster models with {Wo,i^o) = (6, 0.,3) and a power-law mass 
function for power-law index (a) ao = —1.20, (b) qo = —2.35 
and (c) ao = —3.50, respectively. As the slope of mass function 
increases the central velocity dispersion of the most massive stars 
drops faster compared to the other mass function models. The 
behaviour of Uc on pc shows a self-similarity regardless of the 
individual mass. 

5 EVOLUTION AFTER CORE-COLLAPSE 

So far we have concentrated our discussion until the core- 
collapse. We now discuss the evolution beyond the core col- 
lapse. 

5.1 Central density and msiss loss 

Figs. 18 and 19 show the time evolution of the central density 
of the rotating stellar systems beyond the core collapse. The 
evolutions for two component models with the central poten- 
tial Wo = 6 (Fig. 18(a)) and 3 (Fig. 18(b)) are displayed for 
clusters with three different initial rotations {loo = 0.0, 0.3 
and 0.6 for Wo = 6, and uio = 0.0, 0.8 and 1.5 for Wo = 3). 
For two component clusters we employed the mass function 
M2C, while the mass function MCB for the continuous mass 
spectrum is included. The models with the largest initial ro- 
tation (solid lines) reach the collapse earlier than the models 
with smaller initial rotation (dashed and dotted lines). We 
can clearly see the acceleration of the evolution during pre 
and post collapse phases due to initial rotation. 

Fig. 20 displays the run of the individual central den- 
sity for models with the continuous mass spectrum. The 
mass segregation due to the energy equipartition during the 
early stage is shown clearly for all the present models. The 
run of the central density after core bounce roughly fol- 
lows a power- law of pc oc , for all mass groups. The in- 
dex of power-law tends to be smaller for lower mass com- 
ponents. For the non-rotating model we find /3 = —1.54 
for the highest mass group (mio) and /3 = —1.88 for the 
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Figure 15. Evolution of the central angular speed for cluster 
models with {Wo,i^o) = (3,0.8) and a power-law mass function 
for power-law index (a) qo = —1.20, (b) ao = —2.35 and (c) cio = 
—3.50, respectively. Evolutions during the early stage are shown 
since the initial rotation affects largely on the cluster dynamics 
during the early times. 



Figure 17. Evolutions of (a) the maximum rotational veloc- 
ity {Vrot,max) and (b) the maxirrmm specific angular momen- 
tum (Jz.max) at the equator for a model with (Wq, luq, ao) = 
(6, 0.6, —2.35). While Vrot,max of the massive stars increases with 
time (i.e., increasing pc), Jz,max decreases monotonically due to 
the outward transfer of the angular momentum. 




Figure 16. Evolution of the maximum rotational velocity at 
equator for cluster models with {Wo,ujo) = (6,0.3) and a power- 
law mass function for power-law index (a) ckq = —1.20, (b) 
ao = —2.35 and (c) ao = —3.50, respectively. For the steepest 
mass function (ao = —3.50), Vrot.max of the massive stars in- 
crease much than that of the other models with slowly increasing 
mass functions. 



Figure 18. Evolution of the central density for the cluster models 
with the central potential (a) Wo = 6 and (b) Wo = 3 for mass 
function (mg/mi, M1/M2) = (5, 10). The initial number of stars 
in a cluster is iV = 10000 for different initial rotational parameters 
as indicated in the figure. 
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Tidally limited, N= 10000, a„=-2.35, 10-comp. models 
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Figure 19. Same as the Fig. 18 but for models with power-law 

mass function (ao = —2.35). 

lowest mass group (mi). These values are not far from the 
power law slope /3 = —2 of the central density with time 
which was derived by Kim, Lee & Goodman (1998), for 
two-component models. These slopes increase with the ini- 
tial rotation: /3 = —1.83 and —2.42 for the mio and the mi, 
respectively, for model with wo = 0.6. 

We have shown the evolution of the total mass of the 
cluster for all post-collapse models in Fig. 21. The initial 
mass function, the central King's potential and the degree 
of rotation are shown in each panel. The epochs of core- 
collapse are marked by open squares. Note the marked core- 
coUapsc time is obtained not from the present post-coUapso 
models but from the pre-collapse models discussed in §3. Fig. 
21(a) and (b) are the results obtained from two component 
models. The total mass of rotating cluster decrease more 
rapidly than the non-rotating cluster not only pre-collapse 
phase, but also after core-collapse. The raid decrease of total 
mass for rotating model results in smaller tidal radius. 



5.2 Velocity dispersion and angular speed 

The evolution of the mass-weighted central velocity disper- 
sion (cTc) as a function of the central density (pc) is shown in 
Fig. 22. The evolution of the central velocity dispersion af- 
ter core-bounce for all post-collapse models shows a similar 
trend, i.e., the effect of rotation on the central velocity dis- 
persion after core-bounce is very small. Since the core rotates 
like a rigid body, the rotation velocity near the cluster center 
is negligible. While the run of ac on pc during pre-collapse 
shows the power-law behaviour due to self-similarity of col- 
lapsing core for the single mass system (Cohn 1980, Kim et 
al 2002), the stellar system with mass spectrum does not 
show a simple power-law in the early phase as explained in 
Sect. 3.2, due to the mass segregation. The behaviour of ctc 



W„=6, N=10000, a„=- a. 35 
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Figure 20. Time evolution of the central density of the individual 
mass component for models with (Wo, aoi ^) = (6, —2.35, 10000) 
and the initial rotations of (a) cjo = 0.6, (b) luq = 0.3 and (c) ojq = 
0.0, respectively, pc of the massive stars increases very quickly due 
to the development of the mass segregation. 
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Figure 21. Evolution of the total mass retained in a cluster. 
The parameters characterize the cluster model are written on the 
upper-right corner of the each panel. The epochs of core-collapse 
are marked with the open squares. 



© 0000 RAS, MNRAS 000, 000-000 



Dynamical Evolution of Rotating Stellar Systems: III. The Effect of Mass Spectrum 15 




Figure 22. Evolution of the mass-weighted ID central velocity 
dispersion on the central density. The model clusters shown here 
are the same to the models displayed in Fig. 3-21. The behaviour 
of <Tc on pc after core-bounce shows a power-law signature irre- 
spective of the amount of the initial rotation. 



Figure 23. Same as the Fig. 22, but for mass-weighted central 
angular speed on the central density. Only rotating models are 
shown here. A small spike shown in figure (c) is mainly due to 
the numerical error. 



on Pc after core-bounce can be, however approximated as a 
simple power-law even for the multi-mass systems because 
the mass-scgrcgation was already established at the time 
of the beginning of the expansion. The power-law index 7 
can be predicted by applying energy balance argument (Kim 
et al. 2002). Inspection of Fig. 22 shows 7 ~ 0.42, slightly 
smaller than the power-law index of ~ 0.33 obtained for the 
single mass system. 

In Fig. 23, we have shown the run of the mass-weighted 
central angular rotational speed (f2c) on the central density 
(pc) for the whole post-collapse models which have the initial 
rotation. The angular speed also appears to follow roughly 
power law on pc during both pre- and post-collapse, except 
for the early evolutionary stage in the pre-collapse phase for 
the models with Wo = 3. The evolutions of Qc on pc after 
the very early stage shows a similar trend irrespective of the 
degree of the initial rotation. It implies that the effect of 
the initial rotation for the cluster center disappears quickly. 
The evolutions of Uc and f2c on pc of the individual mass 
component for models with the continuous mass spectrum 
is displayed in Fig. 24. The run of CTc is represented well by a 
simple power-law. The higher the stellar mass, the shallower 
the power-law index. While the mass-weighted central angu- 
lar rotation speed appears to show the power-law behaviour 
on Pc Oc of the low mass stars deviates from the power-law 
type evolutionary behaviour during post-collapse phase. It 
shows much steeper decrease than higher mass components, 
though the contribution of the lower mass component on the 
total f2c is negligible. 
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Figure 24. Evolution of the central velocity dispersion (a,b)and 
the central angular speed {c,d) of the individual mass component 
for model clusters with the fastest initial rotation (ojo = 0.6 for 
Wo = 6 and wo = 1-5 for Wo = 3). In ctc on pc plots, the indi- 
vidual mass component show a power-law behaviour, though the 
index of the power-law depends on the mass of the individual star. 
The direction increasing individual mass is denoted with arrows 
in panel (a) and (c), respectively. 
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Figure 25. Radial profiles of the rotational velocities at a 
few selected epochs for a cluster model with (Wo , coo , ao) = 
(6,0.6,-2.35). Half-mass radii (r^) at corresponding epochs are 
marked with open squares. At the initial time, the rotational in- 
side rfi follows a rigid body rotation. Note that the radii is mea- 
sured with the current core radii. 



Figure 26. Same as the Fig. 25 but for a few selected mass 
components, (a) the lowest mass component mi, (b) m4, (c) mr 
and (d) mm, the most massive star. The rotational structure of 
the massive stars extends farther than that of the less massive 
stars. 



5.3 Rotational profiles 

In Figs. 25 and 26, we have shown the radial profiles of the 
rotational velocity at equator for model with (WojCOo) = 
(6, 0.6) and continuous mass spectrum with ao = —2.35. 
The run of the rotational velocity for the selected mass com- 
ponents (mi, m4, my, mio) is displayed in Fig. 26. Epochs 
shown in this plot are t/Trhfi = 0.00,1.50,1.84,3.24 and 
3.90, where the core-collapse occurs at t/Trh,o = 1.84. At 
t/Trh,{) = 3.90 the cluster is nearly completely dissolved. The 
half-mass radii at these each evolutionary times arc marked 
with open squares and the direction with increasing times 
is shown with arrow in Fig. 25. Note that x-axis, i.e., radius 
of cylindrical coordinate is measured in units of the current 
core radius. 

Total rotational velocity decreases continuously with 
time due to the loss of the angular momentum through the 
cluster boundary. The outward transfer of the angular mo- 
mentum causes the increase of the radius when it is mea- 
sured in units of current core radius, where the rotational 
velocity becomes the maximum. The radius of maximum 
rotational velocity decreases after core bounce due to the 
higher mass-loss. While the changes of angular momentum 
due to the outward transfer and the loss through the clus- 
ter boundary changes the structure of the velocity profile 
beyond the half-mass radii, the shape of rotational velocity 
inside the cluster core remains to be rigid body rotation. 



5.4 Evolution of mass function 

Mass segregation drives the concentration of the higher mass 
stars in the cluster center while the lower mass stars go to 
the outer halo of the cluster. This leads to the preferen- 
tial evaporation of lower mass stars through the the tidal 
boundary. Therefore, low mass stars are selectively depleted 
from the initial mass function, resulting a change of the mass 
function with time. Near the end of the evolution, the mass 
function can even be inverted, i.e., there are more high mass 
stars than lower mass stars. The global mass function 0(m) 
at several selected epochs for four different models with the 
continuous mass spectrum is displayed in Fig. 27. We include 
the non-rotating models and the rapidly rotating models in 
order to study the effect of the initial rotation on the mass 
function. The shape of the mass function at the time of 
core-collapse is distinguished by dashed line from the mass 
function at the other epochs. The evolution of power-law in- 
dex {a(M)) for whole post-collapse models with the contin- 
uous mass spectrum is shown is Fig. 28. The mass function 
deviates from the power-law as a result of mass-dependent 
evaporation rates. In Fig. 28, we have shown the typical error 
of the fitted power-law index for the models which has the 
highest initial rotation. The epochs where the core-collapse 
occurs is marked with open squares. In the figures, we have 
chosen M/Mq for the abscissa instead of time because the 
life-time of clusters are significantly different for different 
models. 

It is evident that the behavior of a does not depend 
sensitively on the degree of rotation, although the rotation 
causes less change of a. Although the effect is small, the rotar 
tion causes slower changes in a as a function of M/Mq. The 
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Figure 27. Mass functions of tlic several selected epoclis for the 
cluster models with the continuous initial mass function with a 
power-law index cto = —2.35. For a comparison clusters with- 
out the initial rotation (Figs, a and b) and the clusters with the 
highest rotation (Figs, c and d) for a given central potential are 
displayed. The epochs where the core-collapse happened are dis- 
tinguished with dashed lines from the other evolutionary epochs. 
The change of mass function measured by the power-law index as 
a function of the cluster mass in units of initial mass is shown in 
next figure Fig. 28). 



angular momentum exchange among different mass species 
appears to act against the mass segregation. 

The mass function also becomes dependent on the lo- 
cation as the mass segregation proceeds. The existence of 
initial rotation and its outward transfer also affect the en- 
ergy exchange process that causes mass segregation. In or- 
der to see the role of rotation on the mass segregation, we 
have compared the evolution of a at four different loca- 
tions in Fig. 29. The different locations arc chosen such that 
r < rc, r < Vh, Th < r < 2rh and 2r;j < r < Sr^. The slope 
of mass function inside the core (Fig 29a) increases with 
time during pre-coUapse evolution. At the time of the core- 
collapse the shape of the mass spectrum become inverted for 
all models. The power-law index a at this evolutionary stage 
remains ~ 1.3 with a weak dependence on the the degree of 
the initial rotation. The mass function within the half-mass 
radius is flattened continuously with the decrease of total 
meiss. The model without the initial rotation has a rather 
shallower mass function at fixed total mass compared to that 
of the cluster rotating fast. Initial rotation plays an impor- 
tant role in removing the stars through the tidal boundary 
including both the lower and higher mass stars. Compared 
to non-rotating models, high mass stars in rotating models 
appear to have higher chance of evaporation. Thus the mass 
function changes more slowly for rotating models than non- 
rotating models. Since the effect of the initial rotation re- 
mains in region rn < r < 2rh, it is expected that the shapes 



Figure 28. Evolution of the power-law index Qjvf on the total 
cluster mass for whole post-collapse models with a continuous 
mass function. Open squares represent the core collapse times 
for each cluster model. The errorbar on each panel represent the 
typical errorbar of the fastest rotating model when performing a 
linear least square fitting for a double logarithm plot such as the 
panels in Fig. 27. The slope of the power-law mass function for a 
cluster with the fastest initial rotation is slightly lower than that 
of non-rotating model. 

of mass function in this region shows significant dependence 
on rotation. The concentration of the high mass component 
in the cluster center drives the steepening of the mass func- 
tion near the half-mass radius. When the development of 
mciss segregation is settled in, the slope of the mass func- 
tion decreases with time due to the selective evaporation of 
the lower mass stars throtigh the tidal radius. The evolution 
of a for the highly rotating model shows a different profile 
compared to that of non-rotating model at this radial region. 
While there is a central concentration of the massive stars, 
the angular momentum which has the maximum around rn 
pushes the high mass stars outward. Therefore, there is only 
a little change in a during the early evolutionary stage for 
rapidly rotating models. 



6 CONCLUSION AND DISCUSSION 

We have studied the dynamical evolution of the rotating 
stellar systems with the mass spectrum by solving the orbit- 
averaged 2D FP equation in {E, J^) space. Numerical simu- 
lations are performed both for simple two component clus- 
ters and for clusters with a power-law mass function repre- 
sented by ten mass species. In order to explore the evolu- 
tion after core collapse we add the heating by three-body 
processes. We have employed rotating King models as ini- 
tial models, where the velocity dispersions for all mass com- 
ponents are equal, i.e., no mass segregation at the begin- 



© 0000 RAS, MNRAS 000, 000-000 



18 Eunhyeuk Kim}'^ e-mail: ekim@cfa.harvard.edu, Hyung Mok Lee^ e-mail: hmlee@ astro. snu.ac.kr, and Rainer 




Figure 29. Same as the Fig. 28 but for four different radial re- 
gions of cluster models with the central concentration Wo = 3, 
(a) r <rc, (b) r < r^, , (c) < r < 2/-^ and (d) 2rh < r < 3r,j , 
respectively. A combined effect due to the mass segregation and 
the initial rotation is shown clearly in panels (c) and (d). 



Figure 30. Comparison of observed radial profile of the rota- 
tional velocity over ID velocity dispersion for the galactic glob- 
ular cluster M15 (data provided by Gcbliardt 2002) with the re- 
sult of a selected two-component model with (Wo,uJo) = (6,0.6), 
(m2/mi,Mi/M2) = (5, 10). 



ning. The rotating King models are characterized by two pa- 
rameters: initial central potential (Wo) and degree of initial 
rotation (wo). Clusters with two different central potential 
Wo = 6 and 3 are studied extensively. For models only until 
core-collapse, we have studied the evolution of the clusters 
with six different mass functions for two-component models, 
while three models are studied for clusters with a continuous 
mass function (power-law). For models where the evolution 
beyond the core-collapse are explored, we consider clusters 
with a mass function M2C (m2/mi — 5,Mi/M2 = 10) for 
the two component model and clusters with ao = —2.35 for 
the model with continuous mass spectrum. 

Our results show that, as in equal mass system, the 
presence of the initial rotation accelerates the dynamical 
evolution as manifested by rapid core collapse and dissolu- 
tion. The degree of the acceleration depends on the amount 
of the initial rotation and the shape of mass function. As 
the ratio of the mass (m2/mi) increases the degree of the 
acceleration decreases for two component models since both 
mass segregation and the rotation compete in acceleration of 
core-collapse. If jrvz/ini is very large, the mass segregation 
alone could significantly reduce the core-collapse time and 
there is not much room for rotation to accelerate further. 
For models with the power-law mass function, the acceler- 
ation rate of core-collapse due to the rotation is larger for 
the model with a steeper slope of the mass function for a 
given initial rotation parameter {uiq). The shortening of the 
life-time (dissolution of cluster) due to the rotation is ob- 
served far beyond the core bounce. The increase of mass 
loss rate, resulting from the enhanced two-body relaxation 
process causes the faster dissolution of cluster. 

The evolution of Uc on pc can be approximated by a 



power-law except for the early evolutionary stage regardless 
of the degree of rotation both for pre- and post-collapse. The 
measured power law index 7 for pre-coUapse is very close to 
the value obtained for single mass system, while we have 
obtained a slightly shallower slope for post-collapse phase. 
The development of mass segregation, a consequence of the 
evolution of the multi-mass system is visible clearly both for 
non-rotating and rotating models. The evolution of Q,c on pc 
shows a power-law behaviour, too. The slope of the power- 
law is, however, larger than that obtained for the single mass 
cluster. While the angular momentum is transferred only 
outward for the equal mass system, the exchange of the an- 
gular momentum between different mass species occurs for 
the multi-mass system, resulting a faster increase of fie on 
Pc- Due to a cooperation of the central concentration of the 
massive stars (mass segregation) and the transfer of angu- 
lar momentum from high mass to low mass stars, the radii 
where Vrot reaches the maximum value goes outward. The 
maximum rotation of the most massive stars oven increases 
at the early times, while it shows a monotonic decrease for 
the single mass system. 

There axe a few observations regarding the direct mea- 
surement of the rotation parameter {Vrot) (Gebhardt et al. 
1995). Recently Gerssen et al. (2002) reported the kinemat- 
ical study of the central part of the globular cluster M15, 
including Vrot- M15 is known as the globular cluster which 
contains a collapsed core. The radial profile of Vrot/o (ro- 
tational velocity over one dimensional velocity dispersion) 
shows a steep increase and a rapid decline followed by a 
slow rise. The steep rise of Kot/o" near the cluster center 
can not be explained by a single mass model (Kim et al. 
2002). We have shown the run of Vrot jo for two compo- 
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nent model in Fig. 30 at two selected epochs, one in pre- 
collapse {t/Trhfi = 0.42) and the other after core bounce 
{t/Trhfi = 2.41), together with the observed data by Geb- 
hardt (2002). The radius is measured in current half-mass 
radius. It is evident that cluster loose the angular momen- 
tum through the tidal boundary. The behaviour of the mgiss- 
weighted Vrot/cr (Fig. 30(a)) are roughly similar to that of 
observed profile of M15 except for the central region. Note 
that the known half-mass radius of M15 is 3 .09 (Djor- 
govski 1993). For the high mass component, Vrot,i/oi where 
i represents the individual mass component, aX, t = 2.41 is 
higher than that obtained at i = 0.42, especially beyond th- 
It is mainly due to the lower velocity dispersion after core 
bounce, not due to the higher value of the rotation veloci- 
ties. The highly rotating central region of M15 is, as in sin- 
gle mass system, not explained with the current multi-mass 
models, although there is difference in rotational structure 
between the single mass and the multi-mass systems. From 
the kinematical study of the central region of M15, Gerssen 
et al. (2002) proposed the presence of intermediate mass 
black hole (IMBH) with M = 3.9 x lO^M©. On the con- 
trary, Baumgardt et al. (2002) demonstrated that the rise 
of velocity dispersion into the center can bo explained with 
the clustering of the remnant stars, which are expected to 
be overwhelmingly populated in the cluster core than the 
normal stars. However, they did not rule out the possibility 
of the presence of IMBH. It may be necessary to include the 
remnant stars and the IMBH in current 2DFP models to 
investigate the role of the initial rotation on the observed 
strong increasing of Vrot/cr near the center of M15. 

The models presented here still neglect many important 
physical processes occurring in real star clusters. The stellar 
evolution and primordial binaries arc known to affect the 
early evolution of globular clusters. The possible existence 
of intermediate mass black hole in the center could also affect 
the course of dynamical evolution. These will be the task of 
future works. 

This work was supported by the Korea Research Foun- 
dation Grant No. D00268 in 2001 to HML and by SFB439 
to RS. 
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